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Abstract. In [T] we defined symmetry classes, commutation symmetries and symmetry 
operators in the Hilbert space Ti of the ID spin- 1/2 Heisenberg magnetic ring with N sites 
and investigated them by means of tools from the representation theory of symmetric groups 
Sn such as decompositions of ideals of the group ring C[<Sjv], idempotents of C[5jv], discrete 
Fourier transforms of <Sjv, Littlewood- Richardson products. 

In the present paper we consider symmetry operators a £ C[<Sjv] on Ti which fulfil 

a* — a, (S) 

where a denotes the complex conjugate of a and the element a* £ C[«Sjv] of a = Y2, p a vP 15 
defined by a* = Y2 p a pP • Such symmetry operators are self-adjoint because of the relation 
(au\v) = (u\a*v) which holds true for all a € C[<Sjv] and u, v £ Ti. They yield consequently 
observables of the quantum mechanical Heisenberg model. We prove the following results: 

(i) Let G C Sn be an arbitrary subgroup and x De an irreducible character of G then 
X := ^2 P £G x(p) P i s an idempotent of C[<Sjv] which has property (S). This leads to a big 
manifold of observables. In particular every commutation symmetry belonging to G yields a 
1-dimensional character of G from which we can build one of the idempotents x described here. 

(ii) Let 1Z C C[5jv] be a minimal right ideal of C[5jv]. Then the set of all generating 
(primitive) idempotents of TZ contains one and only one idempotent e which satisfies e* = e. 
This is a result by H. Weyl [2] for which we give a new proof. 

(iii) Every idempotent e with property (S) can be decomposed into primitive idempotents 
e = /i + . . . + fk which have also property (S) and satisfy fi- fj = if i ^ j. We give a computer 
algorithm for the calculation of such decompositions. 

In big group rings C[5jv] computer calculations are only possible by the use of a diskrete 
Fourier transform D : C[Sn] — > © A C dAXdA . We present two algorithms which allow the 
calculation of the matrix B — D\(a*) from a matrix A — D\(a) of an a £ C[Sat] without 
to determine the inverse Fourier transform of A. These algorithms use as precomputed data 
permutation sets V\ which define bases {p • z\\p G V\} of the minimal two-sided ideals Z\ 
belonging to the partitions A h N. {z\ is the central-primitive generating idempotent of Z\). 
The sets V\ can be used for any discrete Fourier transform D. 

In our investigations we use computer calculations by means of the Mathematica packages 
PERMS and HRing. 



1. The Heisenberg model, symmetry operators, symmetry classes, commutation 
symmetries 

We summarize essential concepts of the one-dimensional (ID) spin-1/2 Heisenberg model of a 
magnetic ring (see e.g. (SHU [5] or [IJ). 



We denote by N the set N := {1,... ,N} of the integers 1,2,..., N and by K N the set of all 
functions a : N — > K. 

Definition 1.1 We assign to every function a G K N the sequence \a) := |cr(l), <r(2), . . . ,a(N)) 
of its values over the set N. Then the Hilbert space of a ring model with N nodes and 
a spin-alphabet of K letters {see |6j) is the set of all formal complex linear combinations 

Ti := Cq [\a)\a€ K N | of the \o~) in which the set B := | \o~) \ a G K N X of all \o~) is 
considered a set of linearly independent elements. We equip Ti with the scalar product 

{cr\a') := 8 a ,a' and {u\v) := u a v a for u= u a \o~) , v= v a \a) G H . (1) 

Obviously, Ti is an Hilbert space of dimension dimTi = K , in which B is an orthonormal basis. 
The original ID spin-^ Heisenberg ring arises for K = 2, where o~(k) = 2 represents an up spin 
o~(k) =| and o~(k) = 1 a down spin o~(k) =| at site k. According to [HJ H], the Hamiltonians 
Hp (Ha) of a ID spin-^ Heisenberg ferromagnet (antiferromagnet) of N sites with periodic 
boundary conditions Sf^ +1 := 5f , a G {+, — , z}, are defined by 



N | . 

#4 := -F F , J = const. > (2) 



H F ■= n ( S k S k+l + S k S k+l) + S k S k+l 
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fc=l 

where S 1 ^, S^T and S% are spin flip operators. 

Symmetry classes classes in Ti can be defined in analogy to the definition of symmetry classes 
of tensors (see [21 El El El H] ) . First we define symmetry operators. 

Definition 1.2 Let Ti be the Hilbert space of a ring model of N sites and K letters, Sm be the 
symmetric group of the permutations of 1,2, ... ,N and C[Sn] be the complex group ring of Sn. 

(i) Following \1Q\ 111] we define the action of a permutation p G Sn on a basis vector \o~) G Ti 
by p\a) := \a o p _1 ). 

(ii) Every group ring element a = ^2 p a p p G C[5at] acts as so-called symmetry operator on the 
vectors w = w A a ) ^H by 

aw := ^2 a p w vP\ a ) = a P w <r\ a °P 1 ) ( 3 ) 

A symmetry operator is a linear mapping a : Ti — > Ti which possesses the following properties: 

Proposition 1.3 ([lj) (i) It holds p(q\cr)) = (p o q)\a) for all p,q G Sn, where "o" is defined 
by (poq)(i) :=p(q(i)). 

(ii) It holds a(bw) = (a-b)w for all a,b G C[<Sjv] and w G Ti, where "■ " denotes the multiplication 
a ■ b := ^2 ^2 a p b q p o q of group ring elements. 

(iii) It holds (au\v) = (u\a*v) for all a G C[<Sat] and u,v G Ti, where a denotes the complex 
conjugate of a and the element a* G C[5at] of a = Yl p a pP * s defined by a* = YlpQpP 1 - 

Now we define symmetry classes in Ti. 

Definition 1.4 Let TZ C C[<Sjv] be a right ideal o/C[5tv]. Then 

Tin ■= {au G Ti \ a G TZ , u G Ti} (4) 
is called the symmetry class of'Ti defined by TZ. 



As for tensors one can prove (see e.g. Chap. V, §4] or [12, p. 115]). 

Proposition 1.5 Let e G C[«Sjv] be a generating idempotent of a right ideal 1Z C C[<Sjv], *- e - 
1Z = e • C[Sn] ■ Then a u G TL is in TLn iff eu = u. 

Corollary 1.6 If e, f are generating idempotents of a right ideal TZ C C[5jv]> then u £ TL-r. iff 
eu = fu = u. 

Definition 1.7 We denote by Jq the set Jq := {a £ C[<Sjv] I a« = 0VuG TL}. 

If N > K then Jq ^ {0}, because then the idempotent e := X] p s ig n (p)P satisfies e\a) = 
for all | a) G B. 

Proposition 1.8 ( |12L p. 116]) Jq is a two-sided ideal o/C[<Sjv]- 

Every two-sided ideal of C[«Sjv] has one and only one generating idempotent which is central. 
Let /o be the generating idempotent of Jq. Then / := id — /o is also a central idempotent which 
is orthogonal to fo, i.e. f • fo = fo • f = 0, and which generates a two-sided ideal J := f ■ C[Sn] 
that fulfills C[5/v] = J © Jq. J contains all those symmetry operators a for which ker a C TL. 
These are the symmetry operators in which we are interestet. (Compare [12, p. 116].) 

A special type of symmetries are symmetries of elements u G TL with respect to commutations 
of the nodes of the ring. 

Definition 1.9 Let C C Sn be a subgroup of Sn and e : C — > S 1 be a homomorphism of C 
onto a finite subgroup in the group S 1 = {z G C | \z\ = 1} of complex units. We say that u G TL 
possesses the commutation symmetry (C, e) if cu = e(c)u for all c G C. 

Proposition 1.10 ([12, p. 115]) Let (C, e) be a commutation symmetry. 

(i) Then the group ring element e := Scgc e ( c ) c * s an idempotent of C[C] C C[«Sjv]- 

(ii) A u G TL has the symmetry (C, e) iff e*u = u. 

Consequently, a commutation symmetry (C, e) defines the symmetry class TL-r, of the right ideal 
K = e* ■ C[S N ]. 

2. Self-adjoint symmetry operators 

In the present paper we are interested in symmetry operators a G C[<Sjv] which are self-adjoint 
with respect to the scalar product (•, •) of TL. Such symmetry operators can lead to observables 
of the ring model. First we find 

Proposition 2.1 Two symmetry operators a, b G C[5at] satisfy (au, v) = (u, bv) for all u,v G TL 
iffb = a* + c with c £ Jq. 

Proof. The relations {au,v) = (u,a*v) and (au,v) = (u,bv) for all u,v G TL lead to 
(u, (a* — b)v) = for all u, v G TL. From this we obtain (a* — b)v = for all v G TL. □ 

Lemma 2.2 

(i) If a £ J, then a, a*, a* £ J, too. 

(ii) If a £ Jq, then a, a*, a* £ Jq, too. 

Proof. We show (ii). If au = for all u £ TL, then aw = au = au = for all w = u £ TL, 
i.e. a £ Jq. Further, the condition au = for all u £ TL leads to (au,v) = {u,a*v) = for 
all u,v £ TL. Thus we obtain a*v = for all v £ TL, i.e. a* £ Jq. Finally, a combination of 
a, a* £ Jq yields a* = a* £ J . □ 

Now it follows from Proposition 12.11 and Lemma 12.21 



Proposition 2.3 A symmetry operator a £ J is self-adjoint iff a* = a. 

Proof. According to Proposition 12.11 a is self-adjoint iff a = a* + c with c £ Jq. But since 
a, a* £ J we obtain also c = a — a* £ J such that c = 0. □ 

Definition 2.4 We say that a symmetry operator a £ C[5/v] has property (S) if a satisfies the 
condition 



Every a £ C[5at] with property (S) is self-adjoint. However, an arbitrary symmetry operator 
a = b + b G C[5at] with b £ J ,b £ Jq has only to fulfil F* = 6 to be a self-adjoint operator. We 
investigate in our paper operators a £ C[Sn] which have property (S). 

Now we search for classes of symmetry operators which possess property (S). 

Remark 2.5 Obviously, Young symmetrizers yt £ C[Sn] do not have property (S) in general. 
For instance, one can easily check that the Young symmetrizer yt = [1, 2, 3] + [2, 1, 3] — [3, 1, 2] — 
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[3, 2, 1] £ C[Ss] of the standard tableau t = does not fulfil (S). 

The following theorem yields a big manifold of symmetry operators with property (S). 
Theorem 2.6 Let G C C[5at] be a subgroup and x be an irreducible character of G. 

x ■■= 4rrp ^2x(p)p (5) 
is an idempotent which satisfies (S), too. 

PROOF. The assertion follows from [121 Prop. II. 1.47], = x(p) an d x(id) el. D 

Corollary 2.7 If (G,e) is a commutation symmetry, then e : G — » C is a 1-dimensional 
character of G which is irreducible. Consequently, the idempotent e := r^r ^2 p£ c e (p)P °f eveT "!J 
commutation symmetry has property (S). 

Proof. The assertion follows from [I2j Prop. II. 1.48]. □ 

Already the commutation symmetries lead to a big manifold of symmetry operators with 
property (S). Let us consider the Sq . A list of all commutation symmetries belonging to 
subgroups of Sn with iV < 6 is given in [I2j Appendix A.l]. The Sq possesses 55 conjugacy 
classes of subgroups G / {id} (solvable and non-solvable). The number of solvable subgroups 
G 7^ {id} of ^6 is equal to 1428 (see [T31 III.2.2]). On all theses subgroups one can define 
between 1 and 8 commutation symmetries. This yields a big number of operators with property 
(S). However, if we use irreducible character for the construction of operators with property (S), 
the number of such operators becomes still bigger. 

Theorem 2.8 Let 1Z C C[5jy] be a minimal right ideal with 1Z % Jq. (1Z defines a symmetry 
class). Then the set of (primitive) generating idempotents of 1Z contains one and only one 
idempotent f with property (S). / can be formed from an arbitrary generating idempotent e of 
K by 

f = n e ■ e* , fj, £ R . (6) 

Remark 2.9 A corresponding statement holds true for minimal left ideals C C C[<Sjv]- For left 
ideals one has to replace ([6]) by f = fie* • e. 

1 \G\ denotes the cardinality of G. 



Remark 2.10 A statement similar to Theorem 12.81 was proven by H. Weyl in p. 295]. We 

give here a new proof. 

Proof. Since TZ is minimal and TZ % Jq we have TZ C J". Let e be a (primitive) generating 
idempotent of TZ. Obviously, h := e • e* is a group ring element with property (S). (Note that 
(a • b)* = b* ■ a* .) We show that h is essentially idempotent. 

Since e £ J there exists an u 6 Ti such that eii 7^ 0. Using the scalar product we obtain 
(eu, eu) 7^ and (u, (e* • e)u) / 0. Consequently, we have (e* • e)u 7^ and e* • e / 0. The 
relation e • e* 7^ can be proved when we start a similar consideration from (e*v,e*v), where 
u G W is an element with e*v 7^ 0. 

Now we consider 7^ ((e* • e)u, (e* • e)u) = {eu, (e • e* ■ e)n). This leads to e ■ e* ■ e / 0. If e is 
a primitive idempotent then e • x ■ e is proportional to e for all x £ C[5at]. Since e • e* • e 7^ we 
obtain e • e* ■ e = a e with and h ■ h = e -e* ■ e -e* = a(e ■ e*) = ah. Thus /i ist essentially 

idempotent. 

Next we show that ael. From e • e* • e = a e it follows e* • e • e* = a e*. If we multiply this 
relation from the left by e we obtain a(e* ■ e) = a(e* ■ e) and a = a. So is an idempotent 

which has property (S) because a£l. Finally 1Z' := / • C[«Sjv] = (e • e*) • C[5at] is a right ideal 
with TZ' C 7?.. Since 9 e • e* / and 7£ is minimal we obtain 1Z' = TZ, i.e. / generates 7Z. 

Now we show the uniqueness of /. Assume TZ possesses two generating idempotents /1, 
/2 with property (S). We denote by Tin the symmetry class defined be TZ and by Ti^ the 
orthogonal complement of Tin with respect to (•,•). Every idempotent fi satisfies fix = x 
for all x £ Tin- Since every x S Tin has a structure x = fiZ with z £ H. we obtain 
= (y, fiz) = (fi y, z) = (fiy, z) for all y G Ti^ and all z £ Ti. This leads to fiy = for 
all y £ Tiji- We see that both idempotents f\ and fi have the same effect on Tin and 7^. 
Consequently, we have f\z = fiz for all z £ Ti, from which follows fi — fi £ J§. But since 
TZQj,fi and f 2 lie in J and /1 - / 2 = 0. □ 

3. Decomposition of self-adjoint idempotents 

In this section we present an algorithm for the decomposition of an idempotent e with property 
(S) into pairwise orthogonal, primitive idempotents fi which have property (S), too. 

Theorem 3.1 Every idempotent e £ J with e* = e has a decomposition 

e = fi + f 2 + .-. + fl (7) 

into primitive idempotents fx, . . . fi which fulfill 

f* = fi and fi ■ fj = , i^j. (8) 

Proof. The idempotent e generates a right ideal TZ = e ■ C[Sn] ^ J ■ The decomposition ([7]) 
can be constructed by the following 

Algorithm 3.2 

(i) Form a set y of primitive idempotents y £ C[5at] such that C[5jv] = (B y ey V ' C[«Sjv]- 

(ii) Determine a y\ £ y such that e-y\ 7^ 0. The group ring element e ■ y\ generates a minimal 
right ideal TZ\ := e • y\ ■ C[Sn] with TZ\ C TZ. 

(iii) Determine a primitive generating idempotent h\ ofTZ\ from e ■ y\. 

(iv) Calculate the unique generating idempotent f\ = [i\h\- h\ ofTZ\. 

(v) Calculate the rest r\ := e — f\. 

(vi) Determine further idempotents f 2 , fs, ■ ■ ■ by iteration of the steps (ii), . . . , (v) starting with 
the rests rx,r 2 , ■ ■ ■ instead of e. Stop this iteration when a rest ri = was reached. 



The set y could be constracted from the Young symmetrizers of all standard tableaux of C[<Sjv]- 
In step (ii) the set y is always the original set y. We do not delete from y the elements yi, y 2 , . . . 
found in previous steps (ii). Step (iii) can be carried out by a procedure given in [T2J Prop.I.2.1], 
[T3l Prop.l] or Prop.4.1] 

Now we show that Algorithm 13.21 produces the decomposition ([7]). Obviously, the rest 
r\ = e — /i has property (S). Further it follows e • f\ = f\ from TZ\ C 1Z. This leads to 
ri • /i = 0. Now we can carry out the calculation 

\/x,yeH: = (x,(n- fi)y) = {nxjty} = ((fi-n)x,y) , (9) 

which yield (f\ ■ r±)x = (for all x £ TL) and fx • r\ = since fi,r± E J . From /i • r± = we 
obtain ri • ri = e ■ r\ — f\- r% = r% — = n, i.e. ri is an idempotent with property (S). Thus 
the start of a second iteration of (ii), . . . , (v) with r\ instead of e is correct. The result of this 
iteration are idempotents f 2 , r 2 which fulfil r± = /2 + r 2 , r± ■ fi = fa- 

Consider the symmetry classes H-r.ii H-r^, Tln^ of the right ideals 1Z\ = fi • C[Sn], 
K 2 = h ■ C[S N ], K'i = n ■ C[S N ]. We obtain IZi _L K'i from © and 1Z 2 C K\ from n ■ f 2 = h- 
Consequently, it holds Tt\ _L 7^-2) too. Now a calculation similar to with /2 instead of t\ 
yields f 1 -f 2 = f 2 -fi= 0. 

Since these considerations can be carried out in all iteration steps of Algorithm I3.2I Theorem 
13. II is correct. □ 

4. Use of discrete Fourier transforms 

Computer calculations in big group rings C[5tv] have high costs in calculation time and memory 
because of the high number of permutations in Sjy. One has to use a discrete Fourier transform 
for Sn to reduse these costs. (See [H] for discrete Fourier transforms of groups.) 

Definition 4.1 A discrete Fourier transform for Sn is an isomorphism 



D : C[S N ] -> (g)C dAXdA , D : a= a pP ^ D(a) 

Ah7V pGSjv 

according to Wedderburn's theorem which maps the group ring C[5jv] onto an outer direct 
product &>xi-n C dxXdx of full matrix rings C dxXdx . We denote by D\ the natural projection 
D\ : C[S N ] -> C dxXdx , D x (a) = A x , belonging to the partition N of N . 

The dimension d\ of the d\ x d,\-block matrices ^4^ can De computed from A by means of the 
hook length formula (see e.g. [El QH EH US] or [T2J p.38]). 

The group ring C[5at] decomposes into minimal two-sided ideals Z\, 

c[s N ] = Qz x = Qc[s N ]-z x . (ii) 

AhJV \\-N 

Every Z\ is generated by a unique idempotent z\ which is centrally primitive. z\ can be 
calculated by means of ([5]) from the irreducible character \\ belonging to A h N. Consequently 
every z\ has property (S). 

The structure of block matrices in (I10|) reflects the decomposition (llip since D\ : Z\ — > 
<C dxXdx is an isomorphism. 

Elements a S C[5at] which lie in minimal left, right or two-sided ideals of C[<Sjv], have only one 
non- vanishing block matrix in (jlOf) . In particular this holds true for all primitive idempotents. 
The product a ■ b of such an element a by another element b is in C[5jv] the product of two 



f A Xl \ 

Ax 2 



V 







(10) 



A Xl / 



long sums a = Yl,p a pPi b = YlpbpP of length AM, however in <S) X ^ N C dxXdx only the product of 
two d\ x fi^-matrices A x ■ B x which are much smaller. This leads to a considerable reduction of 
computation costs. 

The steps (i)-(iii) of Algorithm 13.21 were taken from another ideal decomposition algorithm 
which we described in [121 Chap. 1. 2] or [131 19]. The algorithm from |12| , ll3| , l9] runs both in C[<Sjv] 
and in ®ahv < C dxXdx without a necessity of a Fourier transformation between these two spaces 
during the run of the algorithm. The version in <g) XhN C dxXdx is much more efficient than the 
version in C[«Sjv]. 

Because of the relationship between Algorithm 13.21 and the algorithm from [12\ [T5| [§] one can 
transfer Algorithm 13.21 to ®x\-N C dxXdx by means of the remarks in [12, p. 46], [T31 Sec. 5], 
pp.9-10]. 

Only the computation of a* for an element a £ C[«Sjv] can not be transfered to ®ahv C dxXdx 
in a simple way. The calculation of the matrix B = D{a*) from a matrix A = D{a) can not be 
carried out in (£) X ^ N C dxXdx without Fourier transformations between &) X ^ N C dxXdx and <C[«Sjv]- 
In the next sections we present algorithms for the efficient computation of B from A. 

Remark 4.2 For Sn three discrete Fourier transforms are known: 

(a) Young's natural representation of Sn (see [121 P- 51]), 

(b) Young's semi-normal representation of Sn {see [121 p. 55]), 

(c) Young's orthogonal representation of Sn {see [71 pp. 133-135]). 

The DFT (b) is the basis of the fast Fourier transform for Sn by M. Clausen and U. Baum [14j . 
We use (a) as discrete Fourier transform in our Mathematics package PERMS |19| . 
The three DFT map permutations p £ Sn to following types of matrices: 





(a) 


(b) 


(c) 


D{p) 


integer matrix 


rational matrix 


real matrix 



If we want to calculate D{a) of an a £ C[Sn] whose coefficients are only integers or rational 
numbers then (a) and (b) allow calculations not only with a certain numerical precision but with 
infinite precision. In Section [5] we will see that an infinite precision has not an advantage in 
any case. If a calculation lead to such big integers which a computer program can not completely 
decompose into primes than a strong growth of integers can arise and cause memory problems. 

5. Calculation of D(a*) from D(a) 

Now we investigate the question how one can calculate the matrix B = D{a*) G ^Q\\-n C dxXdx 
from the matrix A = D{a) G (^AhN C dxXdx without a discrete Fourier transformation between 
® XhN C dxXdx a ndC[S N ]. 

Proposition 5.1 Every a £ C[Sn] has a unique decomposition 

a = a a where a\ = a ■ z\ £ Z\ . (12) 

Ah 

The a a fulfil 

a* x = a*-z x £Z x . (13) 

Proof. ([12]) follows from (fTTl) . Since z x has property (S) and commutes with every element 
from C[5jv] we have furthermore a* x = z* x ■ a* = a* ■ z x . □ 
Because of Proposition 15.11 we can restrict us to 

Problem 5.2 Determine a matrix-valued function f : <C dxXdx — ► C dxXdx which fulfils 



f(A) =B ^ 3a£Z x : A = D x {a) A B = D x {a*) 



(14) 



Obviously, a function (114ft exists for every Z x . Every such / is a linear function defined on the 
whole of C dxXdx . If we can determine / we obtain a fast algorithm for the computation of B 
from A without a Fourier transformation between ^ X y^ < C dxXdx and C[<Sjv]- 
The simplest algorithm which gives a solution of Problem 15,21 is 

Algorithm 5.3 Let be given a (d x x d\) -matrix A of symbols representing a natural projection 
D\(a) of an arbitrary element a £ Z x . 

(i) Calculate a = (Dx\z x )~ 1 {A). 

(ii) Calculate a* from a. 
(hi) Calculate B = D\(a*). 

B is a (d\ x d x )-matrix of linear functions of the elements of A. The mapping f : A^ B is the 
solution of Problem \h.2\ 

Algorithm 15.31 needs a computer algebra system such as Mathematics which allows symbolic 
computations. It can lead to high costs in calculation time because it carries out a Fourier 
transformations of a relatively small {d\ x cf\)-matrix A to a long group ring element a of length 
N\ and a second Fourier transformation of a group ring element a* of length N\ to a (d\ x d\)- 
matrix B. However these costs arise only once when we run Algorithm 15.31 When we have 
determined the function / then we can very fast compute B from A for all A 6 C dxXdx . 

The costs in memory for Algorithm 15.31 are smaller than for the other algorithms presented 
here. Using Algorithm 15.31 one runs rarely the risk of memory problems. 

A second type of algorithms solving Problem 15.21 bases on the following consideration. For 
every Z x the set {p ■ z x | p £ Sjy} is a generating system of the vector space Z\. We can 
determine a subset V x C <Sjv of permutations such that B\ := {p ■ Z\ \ p G V\} is a basis of Z\. 
If V\ is known we can write every a £ Z\ in the following form 

a = a ■ z\ = a pP " z \ = a P P ' Zx = a ' Zx w ith a := a p p . (15) 

pes N p&Vx ptVx 

It holds true for a, a and z\ 

D x (a) = D x (a), D x (z x )=ld, a* = z* x ■ a* = a* ■ z x . (16) 



Definition 5.4 (i) Let a = J2 p eV x a pP ^ e a 9 rou P r * n 5 element. We denote by v[a p ] the 
column vector in which the a p were arranged according to the lexicographical order of 
permutations. 

(ii) Let A = (A{j) be a matrix. We denote by v[Aij] the column vector in which the elements 
Aij were written line by line. 

Now we formulate our second algorithm which solves Problem 15.21 

Algorithm 5.5 Let A = (Aij) be a symbolic (d x x d x )-matrix and a = Y2peV a pP ^ e a 9 rou P 
ring element with symbolic coefficients a p . 

(i) Determine the coefficient matrix $ of the equation 

v [D x (a)} = $ • v[ap] . (17) 

(ii) Calculate the inverse matrix <3? -1 . 

(hi) Determine the coefficient matrix \& of the equation 

v[D x (a*)] = *-v[a p }. (18) 



(iv) Form the function 

f : viAij] h- • a" 1 • viAij] . (19) 

Arrange the elements of the vectors v[Aij] and • <I>~ 1 -v [Aij] in (d\ x d\) -matrices to obtain 
the function f : C rf * xd * C d * xd \ 

It is clear that Algorithm 15.51 yields the function / searched in Problem 15, 21 Note that 
V\ = dimi^ = d\. So both vector v[a p ] and vector has length d\ and <3? and ^ 

are quadratic (d\ x d|)-matrices. <& has an inverse matrix because B\ is a basis of Z\ and 
-Da|.e a : Z\ — > C dxXdx is an isomorphism. The use of u[o!p] in the right-hand side is correct 
because ct and a* have the same coefficients. Only the places of these coefficients are different 
in a and a*. 

By means of Algorithm 15.51 we calculated the functions / for all partitions A h 5, A h 6 and 
all partitions A h 7, A / (4 2 1), (3 2 l 2 ). The following table shows some information concerning 
the calculations in 57. 



A 


dx 


4 


4 


time 


RAM 


f 


(61) 


6 


36 


1296 


0.6s 




13KB 


(5 2) 


14 


196 


38416 


28.6s 




506KB 


(5 I 2 ) 


15 


225 


50625 


31.2s 




182KB 


(4 3) 


14 


196 


38416 


35.9s 




494KB 


(4 I 3 ) 


20 


400 


160000 


278.8s 


131.7MB 


404KB 


(3 3 1) 


21 


441 


194481 


615.5s 


166.9MB 


2149KB 


(32 2 ) 


21 


441 


194481 


624.3s 


166.9MB 


2187KB 


(3 I 4 ) 


15 


225 


50625 


31.3s 


34.6MB 


182KB 


(2 3 1) 


14 


196 


38416 


36.1s 


34.8MB 


471KB 


(2 2 l 3 ) 


14 


196 


38416 


30.7s 


33.6MB 


433KB 


(21^) 


6 


36 


1296 


0.6s 


14.7MB 


13KB 



The colums of d\ and d\ shows the number of elements in i>[aJ, v [Aij] and ^, respectively, 
'time' and 'RAM' give the computation time and the memory, respectively, which the algorithm 
needed, 'f shows the length of the file in which the function / was stored. 

For the partitions A = (4 21),(3 21 2 ) a big memory problem arose from the use of exact 
integers in the calculation. It turned out, that for instance the matrix $ for A = (4 21) has the 
determinant 

-96974257961965049074024912099372004573795424708044922073363953457157880990336620245014730563217365 
946027469537574499367764896357932972629577610678133530230415476788046519825483493859992633559136 
303784510251434033116082118775724498148053570956977145676159650656795201937451001664676330006575 
437913887075107664190746396865109859747118150041195008145445396420723507449624828927005702309095 
274295858496179822968799982444105484313810897703096007568025010485731184828349379597074706284167 
453662581507561255480875750071503321884946609362947141428351638033559982301650092032000000000000 

000000000000000 
(593 digits) 

During the calculation of <£ -1 rational numbers are generated whose numerators and 
denominators are such large integers as the above determinat. The computer algebra system is 
no longer able to find all prime factors of these integers and to reduce the rational numbers. So 
the integers become longer and longer and the available memory is exceeded. 

We used Algorithm 15.31 for the determination of / for A = (4 2 1), (3 2 l 2 ). The calculation 
had the following characteristics: 



A 


dx 


d 2 
a x 


4 


time 


RAM 


f 


(421) 


35 


1225 


1500625 


4.393h 


440.8MB 


13974KB 


(321 2 ) 


35 


1225 


1500625 









input: Ah iV 
output: V\ 

begin 

i := (xi,x 2 , . . -,x d 2); 

p := id; 

Vx ■= {id}; 

t(0 := H ■ v[D x (p)}; 

solve r(£) =0 =>■ result: x s = cr(rest of the Xj); 
£ := £ <— substitute x s = <r(rest of the Xj); 
while | P A | < d\ do 

p := next-permutation(p); 

r(£) := £ ■ v[D x (p)\; 

if r(£) ^ 

then solve t(£) =0 result: x s = a(rest of the Xj); 
£:=£<— substitute x s = <r(rest of the Xj); 
:=V x U{p}; 

fi; 

od; 

return "Pa; 
end; 



Figure 1. Algorithm for determining of V\ 
6. Determination of the sets V\ 

In this section we present an algorithm which determines the sets V\ for the Algorithm 15.51 (see 
Figure 1). 

We can not simply consider the set {p ■ z\ \ p € 5jv} und determine d\m.Z\ = d\ linearly 
independent vectors in this set, because this set contains iV! vectors from C[<Sjv] of length iV!. 
This data set is very large. 

The algorithm calculates 

D x (p ■ z x ) = D x {p) (20) 

for permutations p £ Sn and checks whether these matrices are linearly independent. 

The algorithm starts with the identity permutation p = id G Sn and V x = {id}. To 
characterize the 1-dimensional linear space spanned by .Da (id), the algorithm uses a vector 
£ = (xi, . . . , x d 2 ) of variables Xj. 

The algorithm forms the scalar product r(£) = £ • v [Da (id)] of the vectors r(£) and [L>A(id)]. 
r(£) is a linear function. The algorithm determines one of the variables x, from the equation 
r(£) = and obtains a relation x s = cr(rest of the Xj). By means of this relation the variable x s 
is eliminated in £. We obtain a vector £ which fulfil £ • w = exactly for the elements w of the 
1-dimensional span of f[-DA(id)]. 

Now the algorithm investigates the next permutation p. To select p, the algorithm uses a 
procedure next-permutation which determines for a given permutation q £ Sn the successor 
p £ Sn according to the lexicographical order of permutations. Then the algorithm checks the 



condition 

r(C) := C ■ v[D x (p)] = . (21) 

If (|2ip is satisfied, then p is canceled. If, however, (|2ip is not fulfilled then a new linear equation 
t(£) = arises which the algorithm uses to eliminate a further variable in £. £ caracterizes now 
the span of two linearly independent vectors u[Z),\(id)], v[D\(p)]. The new found permutation 
p ist added to the set V\. This procedure is iterated until the cardinality of V\ has reached the 
value d\. 

The algorithm needs not much memory because only one matrix D\(p) is stored in the 
computer memory in every step of the algorithm. A halving of the calculation time can be 
obtained by the following theorem. 

Theorem 6.1 Let a set V\ be known for a partition A h N. Then By := {p ■ z x > \ p E Vx} is a 
basis of Zy , where X' h N denotes that partition whose Young frame is transposed to the Young 
frame of A . 

Proof. The generating idempotents z\, zy of Z\, Zy fulfil 

XA(id) v-^ XA'( id ) / \ / ooA 

pG5jv p&Sn 

where xa, Xa' denote the irreducible characters of S n belonging to A, A'. It holds true 
XA(id) = XA'(id), because XA(id) = \/ dim Z\ , xy(id) = V dim Zy and dim^A = dim Zy. 
Further, it follows from [181 Vol. 246, p. 65, Statement 4.11], that 

y p € S N : XX' (p) = sign(p) • x\(p) ■ (23) 

Because of these facts we can write 

z\ = g\ + u\ zy = g\-u\, (24) 

where 

XA(id) v-^ , , XA(id) , . , . 

9\ ■= m l^ Xx ( p ) p > := m\ I^X\{P)P- (25) 



p even p odd 

Now we assume that B\ = {p ■ z\ \ p £ V\} is a basis of Z\, but By := {p ■ zy \ p £ V\} is 
a set of linearly dependent vectors from Zy . Consequently, we can form a non-trivial vanishing 
linear combination of the vectors from By , 

ctpp- zy = . (26) 

P^Vx 

Now the element a := Y^pev x a pP nas a decomposition a = 7 + u, where 7 and u are sums 
which run only through the even or odd permutations of V\, respectively. Now we obtain 

= a ■ zy = (7 + u) ■ (gx - u\) = (7 • - w • u\) + (u ■ g\ - 7 • u\) . (27) 

S v ' V v ' 

= : r =:fi 

It follows from OTD that 



r = , n = o, (28) 

since T and contain only even or odd permutations, respectively. 

Now we consider the non-trivial linear combination (3 := 7 — uj of the permutations from V\. 
The product of (5 and z\ yields 

(3 ■ z x = (7 - u>) ■ (g\ + u\) = (7 • g\ - to ■ u\) + (7 • u x - u ■ g x ) = T - O = . (29) 



So, (3 ■ zx is a non-trivial vanishing linear combination of the vectors from Bx, i.e. Bx would not 
be a basis of Zx in contradiction to the assumptions of Theorem 16.11 □ 



7. An example 

We finish with an example. We consider the Young symmetrizer yt £ C[S§] of the Young tableau 

t= ^ , which belongs to the partition A = (4 1). yt has the length 48. It is proportional 

to a primitive idempotent e = ^ yt and generates a 4-dimensional minimal right ideal of C[5s]. 
Now our algorithms yield 



F :-- 



YS := D x (y t 



YS* := D x {y* t ) 



YS ■ YS* 



1440 



/ \ 



-30 30 

/ 



6 

-6 





1 

■1 




-24 


24 




6 

-6 





_1 1 _1 

4 1 4 

J 



6\ 


-6 

°\ 





D\ was calculated by means of Young's natural representation of . F represents the primitive 
idempotent / with property (S) which generates the same right ideal as yt (see Theorem 12. 8|) . 
/ has the length 120. 

One can easily check the properties YS -YS = SOYS, YS* ■ YS* = 30YS*, F ■ F = F, 
F ■ YS = YS, ^ YS ■ F = F. The relation F* = F can be verified by means of our Mathematica 
package PERMS pi]. 

We implemented all algorithms presented in this paper in PERMS. Mathematica notebooks 
of the calculation for this paper can be downloaded from |20j. 
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